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1 Introduction 



In the chiral limit where the light up, down and strange quark masses go to zero, the QCD lagrangian 
has a SU(3)r x SU(3)l chiral symmetry that is spontaneously broken by the groundstate of the 
theory . There are eight Goldstone bosons: The it's, i^'s and rj. Their interactions are described 
by an effective low energy theory, called chiral perturbation theory (CHPT). At low momenta the 
chiral lagrangian can be expanded in derivatives of the Goldstone fields and in the masses of the 
three light quarks. 

CHPT is a nonrenormalizable theory: Loops produce ultraviolet divergences, which can be ab- 
sorbed by introducing counterterms. The finite parts of the coupling constants of this counterterm 
lagrangian are not determined by the lagrangian that generated the loops. The complete information 
about them is hidden in the QCD lagrangian, but there is no known way to extract this information 
from first principle alone. The couplings have to be evaluated from experiments. As shown by 
Weinberg [Q, the loops are suppressed: Every loop and the associated counterterm correspond to 
successively higher powers of momenta or quark masses. At low energies, the contributions from 
higher loops are small. 

CHPT is an effective theory: It contains all terms allowed by the symmetry of the QCD lagrangian 
in the chiral limit. The local action is generated by the effective lagrangian, which is in general 
not determined by the counterterms alone at higher orders in the expansion of momenta and quark 
masses. Additional terms with finite couplings have to be included. 

The generating functional to one loop and the corresponding local action in the strong sector has 
been determined by Gasser and Leutwyler 0]. An extension to the AS = 1 nonlcptonic weak inter- 
actions to one loop has been given by Kambor, Missimer and Wyler ||. Here we consider virtual 
photons and include them in the mesonic sector. We calculate the divergences of the one-loop func- 
tional in presence of external currents and determine the structure of the local action at this order. 

The article is organized as follows. In the second chapter virtual photons are included in CHPT 
with three flavours. We evaluate the divergences of the generating functional to one loop and de- 
termine the structure of the local action at order p 4 . In the third chapter we discuss the corrections 
to Dashen's theorem at order e 2 m q and compare our estimates with the results in the literature. 
In the fourth chapter we summarize our results and in Appendix A the matrix relations are given 
that are used to simplify the effective lagrangian at order p 4 . Finally, in Appendix B we list the 
renormalized masses of the pseudoscalar mesons at order e 2 m q . 

2 Generating Functional 

We suppose that the reader is familiar with CHPT, otherwise we refer her or him to comprehensive 
reviews {|. 
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At leading order, the effective lagrangian of the strong and electromagnetic interactions which 
respects the chiral symmetry SU(3)r x SU(3)l, P and C invariance is in the mesonic sector [|| 



C 



(Q) 



^F 2 < d^U+d^U + X U + + x + U > + C < QUQU+ >, 



(1) 



where is the field strength tensor of the photon field A^, F^ v — d^A v — d v A^. A is the gauge 
fixing parameter and will from now on be kept at A = 1 . U is an unitary 3 x 3-matrix and incorporates 
the fields of the eight pseudoscalar mesons, 



UU+ = 1, 

U = exp(i$/F ), 

where the A a 's are the Gell-Mann matrices, and 



detU = 1, 
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d^U is a covariant derivative, incorporating the couplings to the photon field A^, the external vector 
and axial vector currents and a M , respectively, 

d M U = d^U - i(v^ + QA^ + a^U + iU + QA^ - a M ), 



where Q is the charge matrix of the three light quarks, 
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(4) 



(5) 



We do not consider singlet vector and axial vector currents and we put therefore tr — tr a M = 0. 
X denotes the coupling of the mesons to the scalar and pseudoscalar currents s and p, respectively, 

X = 2B (s + ip), (6) 

where s incorporates the mass matrix of the quarks, 



M 



m d 



(7) 



F D is the pion decay constant in the chiral limit, F w = F Q [1 + 0{rn q )]. B Q is related to the quark 
condensate < o|mt|o >= — F 2 B Q [1 + 0(m q )], and C determines the purely electromagnetic part of 
the masses of the charged pions and kaons in the chiral limit, 



Ml ± = M 2 K± =2e 2 ^ + 0(m q ). 



(8) 
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To ensure the chiral SU(3) R x SU (3)l symmetry of the lagrangian we introduce local spurions 
Qri Ql instead of the charge matrix Q. The lagrangian is modified, 
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+ -FtKdfU+dpU + X U + +x + U > + C <Q R UQ L U+ >, (9) 
and dfJJ is changed to 

d^U = d^U - i(v^ + QrA/j, + a^U + iU( Vf , + Q L A^ - a M ). (10) 
The rules of chiral transformation are [|[ || 

U g a Ug+, 

Qi -> giQigj, I = R,L 

v^ + QrA^+o,^ -> gR(v f _ l +Q R A^ + a ll )g R ~+ig R d IJ ,g R ~, 

v^ + QlA^- a M -> gLiVft + QLAfi- a^g^ +ig L d fl g£, 

s + ip -> g R (s + ip)g£, 

g R , L G 5(7(3) fl , L . (11) 

Since the currents and a M count as 0(p), where p means the momentum of the external fields, 
the term QA^ has the same dimension. We put 

dim(Q R x) = 0(p), dim(A^) = 0(1), (12) 

the lagrangian £2 ^ s therefore of order p 2 . The convention ([l2|) has the advantage that electro- 
magnetic interactions do not turn upside down the usual chiral counting. The procedure to obtain 
the generating functional at order 0(p i ) is very similar to conventional CHPT. The one-loop graphs 
generated by are of order C>(p 4 ). They contain divergences which are absorbed by adding tree 
graphs, evaluated with the lagrangian C{ Q) of order 0(p 4 ), see below. The generating functional 
becomes, up to and including terms of order 0(p 4 ), 

e i*(« M ,a M) «,p) = N J [ df7 ] [dA ll \e i I* x {^ +c ™} i (13) 

where the integration over the fields is carried out in the one-loop approximation. [Here and below 
we disregard contributions from the anomaly altogether.] To evaluate the divergent part of the 



one-loop functional, we transform the lagrangian (]9[) to Euclidean spacetime, — > C/§ , 

and expand the fields U and A^ around the classical solutions {U^A^) of the equations of motion, 

= U + £-u£u - 2f2< 2 " + ■ ■ ' 
Afi = i M + e M , (14) 
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where we put U — u 2 , and £ is a traceless hermitean matrix. We insert the expansion (|l4|) in the 
action Se, 

S E = [ d 4 X E 



■ J d*x E ji < D^D^- [A„,e][A M ,e] > 



i 

2 

1 ( „ „ I 



-F < (C[£r Kj A^]-H£U M >*m 



^ I + JF? < >J e M | H , (15) 

where the ellipsis denotes higher order terms in the fluctuations £ and e M . We used the expressions 
(see also ref.H) 



T M = hu+ ) d lt u]-liu+G^u-UvG^u+, 



2 L ' M J 2 
^u + d t ,Uu + = --1 



A M = -u + d fJ ,Uu + = --ud^U+u, 



G^ = Vjj, + QlA^ - a M , 

Hij = £/}u + = u + QrU + uQlu + , 

He = u + [Q , U]u + = u + Q nu — uQ eu + , 

(J = -(u + xu + + ux + u). (16) 

r„ and Am are antihermitean matrices, whereas Hrl and a are hermitean ones and / d 4 x E Cf ] 
represents the classical action of cffi . The normal brackets [• • •] denote the commutator, the curly 
brackets {• ■ •} stand for the anti-commutator. We use the parametrization £ = £ a ^ a ; where the 
A a 's are the Gell-Mann matrices. We define a new covariant derivative S M = + X„ with 

= D^ + X^e", 
V = d^ + Xp?, (17) 



where D^ a is understood as 



D»C = d,i a + TfZ b (18) 



and 



rf = -i<[A°,A 6 ]r M >, 
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K P = -X? a = -]f < H L \ a > 5?. (19) 



The action becomes at the one-loop level 



S E \one loop = \J d 4 X E {? {-X^5 ab + ^ + Cl^ 

+£ M (-£„£„ + p) e,,} , (20) 



where now 



<J ab = -\ < [A M ,A a ][A M ,A 6 ] >+\ <a{\ a ,\ b }> 

"iff < [Hr + Hl ' X ^ Hr ~ Hl ' A " ] > "I F ° < HlA<1 ^ >' 
7^ = F < (W A J + X a >, 

P = lK<H 2 L >. (21) 

We collect the fluctuations of the mesons and of the photon field and define a new flavour space r] A , 
where A runs from 1 to 12, rj = . . . , £ 8 , e°, . . . , e 3 ). The action above can now be written as a 
quadratic form, 

S E \one loop = \J d A X E T, A {-^„S AB + A AB ) 7J B , (22) 

where £ M and A AB are 12 x 12-matrices, 

V x; h o 

G ab l^ap 

A = | ). (23) 

This leads to a gaussian integral for the generating functional, 

e -z E \ oneloop = N J dulrfe-i^vv) 

= N'(detV)~i, (24) 
where we defined the differential operator V AB = —T,^T,^5 AB + A AB and 

(/,») = £ jd i x E f A g A . (25) 
Omitting the constant contribution we arrive at 

Z E \ one loop = ^ln(detD). (26) 
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Table 1: P and C transformation properties for the fields, the charges and for their derivatives. The 
defintion of c r Qr and c l Ql is given in the text. Spacetime arguments are suppressed. 



To renormalize the determinant we use dimensional regularization. In Minkowski spacetime, the 
one-loop functional in d = 4 dimensions is given by || 

Z one loop = - ^2^4 / <^ X Tl (t^^^ + \ A J + P 8 ** 8 ' ^ 

where Tr means the trace in the flavour space r\ A and Y av denotes the field strength tensor of , 

Y^ = 0„y„ - d„Y^ + [Y a , Y v ] . (28) 

In table 1 we list P and C transformation properties for the fields and quantities used in the 
one-loop functional. c r Qr and c l Ql mean covariant derivatives of Qr and Ql, respectively, 

c r fl Qi = d fl Q I -i[Gj i ,Q I ], I = R,L, (29) 

that transform under SU(3)r x SU(S)l in the same way as Qr and Ql, 

c lQi -> 9i c'nQi 9t, I = R,L 

g R . L £ SU(3)r. l . (30) 

The divergences in Z on e loop can be absorbed by adding the effective lagrangian that contains 
all terms of order 0(p 4 ) allowed by the chiral symmetry, P and C invariance. Since in the usual 
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physical picture the charge matrix is a constant and has no chirality, we put again Q_r = Ql = Q 
and d^Q = 0. We keep the notation c R ' L Q = —i[G RL ,Q] in order to remember the correct chiral 
transformation of cj^Qj. By using the matrix relations that we list in Appendix A, the lagrangian 
can be simplified to 

Cf ] = L x < d^U+d^U > 2 + L 2 < d^U+d v U >< d^U+dJj > 

+ L 3 < d^U+dJjd v U+d v U > + L 4 < d^U+d^U >< X + U + XU+ > 

+ L 5 < d»U+d^U(x+U + U+x) > + L 6 < X + U + xU + > 2 + L T < X + U - X U+ > 2 

+ L 8 < X + Ux + U + xU + xU + > - iL 9 < d"Ud u U + G R u + d"U + d u UG^ v > 

+ L w < G R ^UG^U+ > + Hx < G R ^G% + G L ^G L ^ > + H 2 < X + X > 

+ KiFl < d^U+d^U ><Q 2 > + K 2 F 2 < d^U+d^U >< QUQU+ > 

+ K 3 Fl (< d^U+QU >< d^U+QU > + < d ti UQU+ >< d^UQU+ >) 

+ K 4 F 2 Q < d^U+QU >< d fi UQU+ > + K 5 F 2 Q < (d^U+d^U + d f "Ud fi U+) Q 2 > 

+ K 6 F 2 < d»U+d^UQU+QU + d>*Ud^U+QUQU+ > 

+ K 7 F 2 a < X + U + U+x ><Q 2 > + K 8 F 2 a < X + U + U+ X >< QUQU+ > 

+ K 9 F 2 Q < (x + U + U+x + XU+ + U X + )Q 2 > 

+ K 10 F 2 o < (x+U + U+x)QU+QU + ( X U+ + U X + )QUQU+ > 

+ K n F 2 < ( X +U - U+ X )QU+QU + ( X U+ - U X + )QUQU+ > 

+ K 12 F 2 a < d^U+ [c R Q, Q]U + d^U [cjQ, Q] U+ > 

+ K 13 F 2 D < c Rf *QUcf;QU+ > + K 14 F 2 Q < c R ^Q c R Q + c L ^Q c L Q > 

+ K lb Ft < QUQU+ > 2 + K 16 F A a < QUQU+ >< Q 2 > + K 17 F A < Q 2 > 2 . (31) 

where G R and G^ u are the field strength tensors of G R , G^, respectively, 

GU = dpGl-hGl-ilG^Gl], I = R,L. (32) 

If we put Q — 0, the terms with the couplings Ki vanish and we are left with the contribution 
to from the purely strong sector, calculated by Gasser and Leutwyler jl. In ([H]) we have 

only considered contributions made out of the building blocks in table 1 . In particular, terms 
proportional to polynomials in are omitted. The coupling constants Li, Hi and Ki are defined 
as 

U = TiX + LKf,), 
Hi = AiX + HKfj,), 

Ki = ViX + KKjj). (33) 
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A contains a pole in d — 4 dimensions, 

A^I^-i^ + nD + l]}. (34) 

The coefficients I\, Aj and the finite parts L^fx) are listed in ref.Q, whereas H[ and iJ£ are of no 
physical significance. The coefficients Ej are 



sx = ! 


E 2 = Z 




E 4 = 2Z 


E 5 - -4 


E 6 = \Z 


E 7 = 


E 8 = Z 


= — 4 


Eio = 4 + \Z 


En = \ 


S12 = 4 


E 13 =0 


E 14 = 


E15 = | - 



(35) 



3Z + 20Z 2 



Ei 6 = -3- IZ-AZ 1 



E 17 = §- §Z + 2Z 2 



with 

^ = ( 36 ) 

r o 

We choose the renormalized coefficients 1£J (/1) in such a way that the Ki themselves do not depend 
on the scale /i, 



/i^ = E i ^-J+ M A^ (/() j, = „. ,:57) 



3 Corrections to Dashen's Theorem 

This chapter is separated in four sections. In the first section we calculate the formal expression 
for the corrections to Dashen's theorem at order e 2 m q , i.e. the squared mass difference (M^ (± — 
M^a) — (Af 2 ± — M^o). In the second section a numerical estimate is given with an upper limit on 
the unknown terms. In the third section we give an independent approach on the basis of ref. 0. 
Finally, we compare in the fourth section the obtained values with the results in the literature. 

3.1 Masses at Order e 2 m q 

In the chiral limit, the masses of the pions and the kaons are of purely electromagnetic nature. A 
relation between the squared masses at order e 2 is given by Dashen's theorem j8|, 

(Ml ± - Ml a ) e m = (M|. ± - M| ) e m (38) 

with 

M 2 ± = A4 ± =2e 2 -|, 

O 

M 2 = M% = for m q = 0. (39) 
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At lowest order in the quark mass expansion the squared masses of the pseudoscalar mesons are 
denoted by M|>. We neglect the mass difference — m u and replace m u , by to — (to u + m^) /2. 
If we expand the lagrangian in the parametrization (|2|) and use the representation (^) of the 
physical fields, the squared masses are 





2e 2 ^ 


+ 2toB q , 


Ml 


= Ml = 


2toB g , 


M 2 K± 


„ 2 c 

= 2e V 


+ (to + to s )-B, 




= M% = 


(to + m s )B , 


ill 


= -(m4 


- 2m s )B a . 



(40) 

This decomposition depends on the convention adopted for the quark self-energies. The effect of 
this ambiguity is small and neglected in this article. The correction to Dashen's theorem at order 
e 2 m q is denoted by 

AM 2 , - AM 2 = (M% ± - M^o) - {M 2 ± - M 2 a ) for m u = m d =m. (41) 

We calculate the Fourier transform of the two-point function of the mesons to one loop, keeping the 
fields associated with the Gell-Mann matrices A Q instead of the physical ones. For the charged fields 
it has a cut at p 2 = M 2 , 

tx e*" < o\TM^a(o)efM c< ^(y)+cfHy)} lQ > ^ ^ 

Za(Mi) 



M a 2 (l -p 2 /^a) 1+e2/o(p2) 



a = 1,2,4,5 (42) 



with 



whereas for the neutral fields it contains a pole, 

d 4 x <o|T^(x)^(o) e l /^{ £ -- (! ' )+ ^ Q,to) }|o> | toone;oop 

Za(M 2 ) , 



Ml -p 2 



a = 3, 6, 7, 8. (44) 



C-2int represents the interaction part of the lagrangian Z a (M 2 ) and Z a (M 2 ) are due to the 

field renormalization and the ellipses denote terms which are not relevant for the mass-shifts at order 
e 2 m q . The contributions to the two-point function are shown in fig. 1 : (b) and (c) are one-photon 
loops. The first contributes to the two-point function of the charged fields only, the latter vanishes 
in dimensional regularization. (d) is the tadpole with a four-meson coupling and (e) represents the 
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contributions from the lagrangian £4 . The results for the different mesons are explicitly given in 
Appendix B. Here we consider only the difference AM|- — AM^, 



AM# - AM; = -e 1 



16tt 2 



3M 2 In ^ - 3M 2 In ^ - 4(M| - M 2 ) 



C 



-16e 2 ^( M )(Af 



A/f 2 A/f 2 
M £ln^-M 2 (21n^- 

/IT /i 2 



Mi 



1) 



+0(e 2 m 2 ), 



(45) 



where L^(p) and -R^k are the contributions from £ 4 

2 



(Q) 



3iTJ + 2K 



Kl - 6-fcT 



10^0 " 12^n) 



(46) 



At the end we need L£ and seven coupling constants K.\ from to determine the correction to 
Dashen's theorem at order e 2 m a . 



3.2 Numerical Results 

We put F equal to the physical pion decay constant, F^ = 93.3 MeV and the masses of the mesons 
to Mtt = 135 MeV, Mk = 495 MeV. C can be expressed as an integral over the difference of the 
vector and axial vector spectral functions as established by Das, Guralnik, Mathur, Low and Young 
Q . If we consider the low-lying vector and axial- vector mesons p and A\ , respectively, and use the 
Weinberg sum rules to eliminate the parameters of the A\, we arrive at jl]J 



C 



^M 2 F 2 ln 



F 2 

p 



p2 _ p2 



(47) 



where M p is the mass of the p meson, M p = 770 MeV, and F p denotes the p decay constant, 
Therefore 



F p = 154 MeV 



C = 61.1 x 10~ 6 (GeV) 4 



(48) 



At tree level, the difference of the squared pion masses is 

AM 2 



(M 3r ±+M 7r o)(M^±-M 7r o) 
2M V x 4.8 MeV, 



(49) 



which is very close to the experimental value, (M f ± — M n o) exPi = 4.6 MeV [Q. The spectral 
functions can also be extracted from data of the r-decay [|l3|. However, this method contains 
uncertainties because the spectral functions can only be evaluated for momenta p 2 < 2(GeV) 2 and 
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therefore we keep our approach above. L^(n) is taken from table 1 of ref.||. Its values at the scale 
points (i = (0.5 GeV; M p ; 1 GeV) are 

LlQi) = (2.4; 1.4; 0.8) ± 0.5 x 10~ 3 . (50) 

These values lead to the following result for the squared mass differences, using again the three 
scales fx = (0.5 GeV; M p ; lGeV), 

AM% - AM 2 = (-0.26; 0.52; 0.99) ± 0.13 x 10~ 3 (GeV) 2 + R^M 2 + R K M 2 K . (51) 

There is a strong dependence on the scale /z for all the terms. Of course, this scale-dependence 
cancels in the full expression. 

Next we calculate the ratio (AA/^/AAf 2 ). For this purpose we put AAf 2 = (M 2 ± - M% ) exp . and 
arrive at 

AMI , . R^MI + RkMI , x 

J = (0.79; 1.42; 1.80) ± 0.11 + ^ • (52) 

In we can estimate the coupling constant K%{p) if we neglect the unknown contributions from 

£4^ proportional to AT 2 (see Appendix B), 



a 1 C 
" 6 8^ Fl 



-16e 2 -| [(Af 2 + 2M» )LJ + M»Lg] + 8e 2 M 2 K K r s + 0(e 2 M 2 ). (53) 

" o 

The values of £4^) are at the three scalepoints fi = (0.5 GeV; Af p ; 1 GeV), 

LXji) = (0.1; -0.3; -0.5) ± 0.5 x 10~ 3 . (54) 
We put AM 2 equal to the observed value and Kg(fi) becomes 

k E(m) = 4 -°; 5 - 6 ) ± L7 x 10-3 ■ ( 55 ) 

Kg is of the same order as the coefficients L\ Q and as it was expected by the so-called naive chiral 
power counting [jl4|. Note that K^(p) has not the usual scale-dependence implied by (j37|), because 
we have dropped scale-dependent terms in (53). Now we assume that the -KJ's which contribute to 
Rn and Rk have the upper limit 

1**1 ~ TT~2 = 63 x 10_3 • ( 56 ) 

lD7T Z 

The contributions from R^M 2 + R^M\ to the squared mass difference AM\ — AM 2 is smaller 
than 

\R V M 2 + R K M 2 K \ ~ 2.6 x 10" 3 (GeV) 2 , (57) 
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which is a large number. Similarly the contribution to the mass ratio (AM^/AM 2 ) is 

\R*M* + R K M$\ < 2 



AM 2 

On the basis of (|56"|), large corrections to Dashen's theorem can therefore not be excluded. 
3.3 Independent Approach 

We give a further estimate by looking at the QCD mass difference of the kaon. In this calculation 
we put irid — rn u ^ 0, but neglect corrections of order e 2 (m,d — m u ). The observed squared mass 
difference of the kaon can be divided into two parts, 

(M 2 K± - M 2 KO ) exp . = (M 2 K± - M 2 Ko ) QCD + (M 2 K± - M% a ) e . m , . (59) 

The same is true for the pions, but the QCD term is of 0[(md — ra u ) 2 ], which is very small, 
numerically (M„± - M^)qcd = 0.17MeV §]. We put therefore 

(M 2 ± - M 2 a ) exp . = (M 2 ± - M 2 ) e . m . . (60) 

The electromagnetic part of the squared mass difference using the notation of ([ll]) becomes 

AM 2 K ~ AM 2 = (AM 2 K ~ AM 2 ) exp . - (AM 2 K ) QCD + 0{e 2 m 2 q ). (61) 

The experimental part is (AM| - AM 2 ) exp . = -5.25 x 10- 3 (GeV) 2 @. We refer to Leutwyler § 
for the calculation of (AM^)qcd without using Dashen's theorem, 



(AM 2 K ) QCD = -{M 2 K - M 2 S^-^ [1 + A M + 0(m 2 )] . (62) 

TO, — TO 



R = (m s — m)/(md — m u ) is the ratio determined from the mass splittings in the baryon octet and 
from p — lo mixing JTT| with the value R = 43.7 ± 2.7. One way to obtain the correction Am is 
an estimate of r\ — rj mixing. The phenomenological informations [l5| indicate that the mixing 
angle is somewhere between 20 and 25 degrees. The range \Q m ' \ = 22° ± 4° corresponds to a value 
of Am = 0.0 ± 0.12. Inserting these uncertainties into the equation above, the QCD squared mass 
difference of the kaon becomes 

(AM 2 K ) QCD = -(5.19 ± 1.01) x 10- 3 (GeV) 2 , (63) 

and the electromagnetic part of the squared mass difference is therefore (we neglect the errors in 
the experimental values) 

AMf - AM 2 = -(0.06 ± 1.01) x 10" 3 (GeV) 2 (64) 
with a large uncertainty for a small value. Similarly the mass ratio (AMS-/ AM 2 ) is 

— — | = 0.95 ± 0.81 . (65) 
AM 2 v ; 



12 



3.4 Comparison with Other Results 

Maltman and Kotchan Jl6| used the effective lagrangian and calculated the one-loop contribu- 
tions to the squared mass differences AM 2 K . They consider terms proportional to e 2 M 2 ln(M 2 /pi 2 ) 
and e 2 Mj c In(M|-/// 2 ) at the scale point /U = 1.1 GeV. The scale independent terms e 2 M 2 and e 2 M] ( 
as well as the contributions from the effective lagrangian of order p 4 are neglected. Their result is 

(AM 2 , - AM 2 ) log = 0.76 x lCT 3 (GeV) 2 , (66) 
which agrees with the contribution from the logarithms in our calculation at [i = 1.1 GeV, 

{AM 2 K - AM 2 )i og = 0.77 x 10~ 3 (GeV) 2 . (67) 
The authors expect that the ratio (AM^ / AM 2 ) is in the range (at this scale point) 

— — | = 1.44 ±0.20. 68 

AM 2 y ' 

Using again AM 2 — (AM 2 ) exp ,, this corresponds to 

AM 2 K - AMI = (0-55 ± 0.25) x 10~ 3 (GeV) 2 . (69) 

Recently two other estimates have been made by Donoghue, Holstein and Wyler jl7) and Bijnens 
fist . The first work is based on chiral symmetry and vector meson dominance. Their result shows 
a total ratio which is rather large, 

- 1-8 ■ (70) 
AM 2 v ; 

The value for AM 2 they obtain at this order is about 20% higher than the experimental value, 

M 2 ± - Af 2 = 2M W x 5.6 MeV. (71) 

Inserting ( |7l|) into the ratio above (|70|), the electromagnetic part of the squared mass difference 
becomes 

AM 2 K - AMI = L23 x 10~ 3 (GeV) 2 . (72) 
This result is in good correspondence with the one obtained by Bijnens p8| , 

AM 2 - - AM 2 = (1.3 ± 0.4) x 10" 3 (GeV) 2 , (73) 

who calculated this value by splitting the contributions to the masses into two parts. A long distance 
part is estimated using the (l/iV c )-approach and a short distance contribution is determined in terms 
of the coupling constants L\ of the lagrangian . 

The results we listed in this section show rather large corrections to Dashen's theorem. They coincide 
only partly with the result we obtained in (|64|), but they could be in agreement with the estimate 
on the basis of (|56|). In order to establish them, we would have to numerically evaluate the relevant 
coefficients K\ of . 
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4 Summary 

1. We have evaluated the divergent part of the generating functional of CHPT including virtual 
photons to one loop. We have calculated the structure of the local action at order p 4 with 29 
coupling constants {L\, . . . , Lio; Hi, Hz; Ki, . . . , Kyj), which have in general a single pole in d — 4 
dimensions. 

2. In a next step we have calculated the one-loop contribution to the meson masses in the isospin 
limit m u — m,d — rh and extracted the correction to Dashen's theorem lp[, 



AMI - AM? = -e 2 



1 



2 1 C 

C 



3M| In - 3Af 2 In ^ - 4(M 2 - M 2 ) 
M|ln^f _M 2 (21n^ + l) 



-l&e 2 —U 5 ^){M 2 K - Ml) + i?.(,i)Af 2 + R K {^)M 2 K 

" O 

+0(e 2 m 2 ), (74) 

where AMp = Mp ± — M po with P the pseudoscalar meson in question. L r 5 (y) and R v ,k are the 
contributions from 

Rn = \e 2 (6K^ - 3KI + 2KI - \W r w - l2K r n ) , 

Rk = -^{Ki + Kl-Mlv-Wl,). (75) 

At the three scale points /i = (0.5; 0.77; 1) GeV we got numerically 

AM| - AM 2 = (-0.26; 0.49; 0.99) ± 0.13 x 10" 3 (GeV) 2 + i? w M 2 + R K M 2 K . (76) 

We assumed the upper limit of R n M% + RkM\ to be 

\RkM 2 + R K M 2 K \ ~ 2.6 x 10~ 3 (GeV) 2 , (77) 

from which we can not exclude large corrections to Dashen's theorem. Note that in the numerical 
part in (|7^) the correction due to £5 (a*) is already included. 

3. In addition, we have given an estimate of AM| — AM 2 on the basis of ref. 

AMf - AAf 2 = -(0.06 ± 1.01) x 10" 3 (GeV) 2 , (78) 
a small number with a large uncertainty. 

4. We have compared our values to the results in the literature [|l6| [l7], We found that 
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they coincide only partly with the result (|7q), but they could be well in agreement with the esti- 
mate in ( |76| ) . In order to complete the expression for AM| — AM, we need the numerical evaluation 
of Rtt^k, i-e the calculation of the finite parts of the coupling constants Ki. This is beyond the scope 
of this work. 
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A Matrix Relations 

For the derivation of the effective lagrangian C± at order p 4 , we used the equation of motion 
obeyed by U, 

d^d^U+U - U+d^dfJJ + (U+x - X + U) 

+^(U + QUQ - QU+QU) - i < U+x - X + U >= 0. (79) 

The Cayley - Hamilton theorem for a 3 x 3-matrix A, 

A 3 - < A > A 2 + - {< A > 2 - < A 2 >} A 

-- {< A > 3 -3 < A >< A 2 > +2 < A 3 >} 1 = 0, (80) 

implies a trace identity for 3 x 3-matrices A, B, C and D, 

< ABCD > = - < ABDC > - < ACBD > - < ACDB > 

- < ADBC > - < ADCB > + < AB >< CD > 

+ < AC >< BD > + < AD >< BC > + < A >< BCD > 
+ < A >< BDC > + <£>< ACD > + <B>< ADC > 
+ <C>< ABD > + <C>< ADB > + <D>< ABC > 
+ <£>>< ACB > - < A>< B ><CD > 

- < A>< C >< BD > - < Ax D >< BC > 

- < B >< C >< AD > - < B >< D >< AC > 

- <C >< D >< AB > + < Ax B ><C >< D > . (81) 

Furthermore we used the observation that Q 2 can be written as 

° 2 = i Q+ ¥ 1 - (82) 
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We confirmed the 17 independent terms < • • • > (i = 1,2,..., 17) in calculating the 

contributions from tree graphs at order p 4 to 

• the constant term e 4 l, 

• the masses M 2 ± ,M 2 and M 2 , 

• the scattering amplitudes at order e 2 m q 

tt + it~ — ► 7T + 7T~, 

-» if ^ , (83) 



the scattering amplitudes at order e 4 



7r + 7r — > ir + n , 



7T + 7T 



if°i4r°, (84) 



the matrix elements 



< 7r + |u75C?|7r + 7r >, < ir + \11j5d\r1r1 >, 

< 0|u7 M 7 5 d|7r _ >, < 0\T u(x)^j 5 d(x)d{y)j^ 5 u{y)\0 >, 

< 7r+\Tu(x)^d(x)d(y)^u(y)\n+ >, < 0\^x)^d{x) d(y) lt ,u(y)\0 

B The Formal Expressions for the Masses at Order e 2 m q 

The photon loop (see fig.l) gives a contribution to the masses of the charged particles, 

1 / M 2 \ 

5M 2 c± \^ loop = -e 2 —M 2 (31n-f - 4) , (86) 

where C stands for 7r and K, respectively. 

From the tadpole, the masses of the charged particles get a change of the form 



M 2 M 2 
A c± M 2 \n^ + B c± M 2 K \n^f- 



(87) 

with (An±, B 7T ±) = (2, 1) and (A K ±, B K ±) = (1, 2). The contributions to the masses of tt° and rj 



are 

i c r / M 2 \ / M 2 \ ' 

SM 2 \ tadpole = e 2 ^^ \a P Ml (in + lj + (3pM% (j P In -f + lj j , (88) 

with the coefficients (a^, /?„•<>, 7^0 ) = (6, 0, 0) and (a,,, /?,,, 7,,) = (—2, 1, 4). The masses of the 
neutral kaons do not get contributions from the loops. At order e 2 m qi the mass-shift due to 
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involves in general only terms with the couplings Ki, but for the charged particles, where L4 and 
L 5 enter as well, 

5M 2 c± \ ciQ) = -16e 2 ^ [(Ml + 2M 2 K )L\(p) + M 2 L r 5 (^)] +■■■. (89) 

4 " o 

The term proportional to L\(p) drops out in the difference AM|- — AM 2 . Explicitly, the masses 
are at 0(e 2 m q ), 
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Ml ± _ e 2 _J_ Ml ( 3 i n Ml _ 4 
16tt 2 x V ^ 2 



|2 



2M 2 ln^+M 2 r ln^ K 



M 2 i 

2 



2 1 C 

~ 6 8^2 

-16 e 2 ^[(M2 + 2M^ + M^] 

^ o 

~e 2 Ml (6X[ + &K r 2 + 5XJ + 5X 6 r - 6X£ 



-15ifJ - 5X 9 r - 23X[ - 18^) 



+8e 2 M 2 K K r 8 , 



-e 2 M 2 (12X[ + 12X 2 r - 18X5 + 9X£ + 10X 5 r 



+iox<r - i2X£ - uk; - iok; - 10K[ { 



10) > 



M| ± -e 2 T ^M|(31n^-4 



— e 



1 C 

8tt 2 

C_ 

p4 



M 2 M 2 
M 2 ln^ + 2M|ln^f 



_l 6e 2_ [ (M 2 + 2A ^ )L r + M | L g 

+ ^e 2 M 2 (3X 8 r + X 9 r + X[ ) 

~^e 2 M 2 K (6X1 + 6X 2 ' + 5X£ + 5X 6 r - 6Xf 



-24X£ - 2X£ - 20X[ - l&K^) 



M 



K" 



Mh ~ (3X[ + 3X 2 + Kl + XJ 

-3xj - mi - Kl - X[ ) : 



M 2 - e 2 — !- — 
" 48^ 2 



2Mf In 



Mi 



| 41n- 



M 



K 



+ 1 
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+ ^Ml{Kl + Kl Q ) 

~e 2 M 2 {12K{ + 12K r 2 - QK r s + 3KI + 6K; 

+6K£ - l2K r 7 - \2K r s - 4KJ - AK r w ) . 



The squared mass differences of the pions and the kaons are 



AM; = 



M 2 ± M 2 a 

2 4- e2 i^ 2 ( 3ln f - 4 ) 



— e 



1 C 
8tt 2 
C 



> 1 -'-J K iiJ 



-16 e 2 - ¥ [(M2 + 2M^)^ + M^] 

^ o 

+2e 2 Af 2 (-2^ + if 4 r + 2X 8 r + 4tf[ + AK^) 
+8e 2 M 2 K K r 8 , 



A' 



M 2 K± - M 2 



A"" 



2e 2 ^ - e 



F 2 



M 2 



T ^M|(31n^f-4 



,2 



M 2 In + 2Mx In 



2 1 C 

-16e 2 ^[(M 2 + 2M|)^ + M|^] 

^ o 

+ ^ e 2 M 2 (3^ + ^ + X[ ) 

-^M 2 - + ^ - QKl - 6K r w - 6K r u ) . 



And finally the difference AM\ - AM 2 is 

o 1 



AM| - AM 2 



3M 2 K In ^ - 3M 2 In ^ - 4(M| - M 2 ) 

A/ 2 A/ 2 
M|.ln^f -M 2 (21n% + 1) 



16tt 2 

2 1 C 

6 8tt 2 F4 

16e 2 ^^( M )(A4 - M 2 ) + R^)M 2 + R K (^)M 2 K , 



where R„ t K are, 



R n = -e 2 (6K r 3 - 3K% + 2KI - 10if[ - \2K\ X ) 
Rk = -U 2 (K r 5 +K r 6 -6K[ -6K r n ). 
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